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Specification 
Title of the Invention 
Electromagnetic Flowmeter 

5 Background of the Invention 

The present invention relates to an 
electromagnetic flowmeter which measures the flow rate 
of a fluid to be measured, which flows through a 
measuring pipe and, more particularly, to an exciting 

10 method and signal processing method capable of realizing 
accurate flow rate measurement . 

An electromagnetic flowmeter measures the flow 
rate of a conductive fluid to be measured, which flows 
through a measuring pipe, by converting the flow rate 

15 into an electrical signal by using electromagnetic 
induction. Fig, 11 shows the arrangement of a 
conventional electromagnetic flowmeter. 

This electromagnetic flowmeter comprises a 
measuring pipe 11 through which a fluid to be measured 

20 flows and a pair of electrodes 12a and 12b which oppose 
each other in the measuring pipe 11 to be perpendicular 
to both the magnetic field applied to the fluid to be 
measured and an axis PAX of the measuring pipe 11 and 
also come into contact with the fluid to be measured. 

25 The electrodes 12a and 12b detect an electromotive force 
generated by the magnetic field and the flow of the 
fluid to be measured. 
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The electromagnetic flowmeter also comprises 
an exciting coil 13 which applies, to the fluid to be 
measured, a magnetic field perpendicular to both the 
measuring pipe axis PAX and an electrode axis EAX that 
5 connects the electrodes 12a and 12b, a signal conversion 
section 15 which detects the electromotive force between 
the electrodes 12a and 12b, and a flow rate output 
section 16 which calculates the flow rate of the fluid 
to be measured from the interelectrode electromotive 

10 force detected by the signal conversion section 15. 

In the electromagnetic flowmeter shown in 
Fig. 11, a plane PLN which includes the electrodes 12a 
and 12b and is perpendicular to the direction of the 
measuring pipe axis PAX is defined as a boundary in the 

15 measuring pipe 11. At this time, symmetrical magnetic 
fields are applied to the fluid to be measured on both 
sides of the plane PLN, i.e., the boundary in the 
measuring pipe 11. The exciting coil 13 can be excited 
by a sine wave exciting method or a rectangular wave 

20 exciting method (e.g., "A to Z of Flow Rate Measurement 
for Instrumentation Engineers" edited by Japan Measuring 
Instruments Federation, Kogyogi justusha, 1995, 
pp. 143-160). 

The sine wave exciting method that uses a sine 

2 5 wave as an exciting current for an exciting coil is 
readily affected by commercial frequency noise. 
However, this problem can be solved by a high-frequency 
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exciting method which uses an exciting current having a 
higher frequency. 

The high-frequency exciting method is 
resistant to 1/f noise such as electrochemical noise or 
5 spike noise. In addition, this method can improve the 
response (a characteristic which makes a flow rate 
signal quickly follow a change in flow rate). 

However, the conventional sine wave exciting 
method is readily affected by in-phase component noise. 
10 An example of in-phase component noise is a shift of the 
amplitude of a magnetic field applied to a fluid to be 
measured. 

In the conventional electromagnetic flowmeter, 
when the amplitude of the exciting current supplied to 

15 the exciting coil varies (shifts) due to a fluctuation 
in power supply voltage, and the amplitude of the 
magnetic field applied to the fluid to be measured 
shifts, the amplitude of the interelectrode 
electromotive force changes, resulting in a flow rate 

20 measurement error due to the influence of shift. Such 
in-phase component noise cannot be removed even by the 
high-frequency exciting method. 

To the contrary, the rectangular wave 
exciting method that uses a rectangular wave as an 

25 exciting current to be supplied to an exciting coil is 
resistant to in-phase component noise. 

In the rectangular wave exciting method, 



- 3 - 



however, the interelectrode electromotive force is 
detected when a change in magnetic field becomes zero. 
When the exciting current has a high frequency, the 
detector must have high performance. 
5 Additionally, in the rectangular wave exciting 

method, when the exciting current has a high frequency, 
effects of the impedance of the exciting coil, the 
exciting current response, the magnetic field response, 
and an overcurrent loss in the core of the exciting coil 

10 or measuring pipe cannot be neglected. It is difficult 
to maintain rectangular wave excitation. 

As a result, in the rectangular wave exciting 
method, high-frequency excitation is difficult, and an 
increase in response to a change in flow rate or removal 

15 of 1/f noise cannot be realized. 
Summary of the Invention 

It is an object of the present invention to 
provide an electromagnetic flowmeter which can remove 
in-phase component noise and correct any flow rate 

20 measurement error and also realize high-frequency 
excitation . 

In order to achieve the above object, 
according to the present invention, there is provided an 
electromagnetic flowmeter comprising a measuring pipe 

25 through which a fluid to be measured flows, an electrode 
which is arranged in the measuring pipe and detects an 
electromotive force generated by a magnetic field 



- 4 - 



applied to the fluid and flow of the fluid, a first 
exciting coil which is arranged separately from a plane, 
which includes the electrode and is perpendicular to a 
direction of an axis of the measuring pipe, and applies 
5 a first magnetic field having a first frequency to the 
fluid, a second exciting coil which is arranged on a 
side opposite to the first exciting coil with respect to 
the plane and applies, to the fluid, a second magnetic 
field obtained by amplitude-modulating a carrier having 

10 the first frequency by a modulated wave having a second 
frequency, a power supply section which supplies an 
exciting current to the first exciting coil and the 
second exciting coil, a signal conversion section which 
separates a component of the first frequency from the 

15 electromotive force detected by the electrode to obtain 
an amplitude, separates one of components of sum and 
difference frequencies of the first and second 
frequencies from the electromotive force to obtain an 
amplitude, and obtains a ratio of the amplitudes, and a 

20 flow rate output section which calculates a flow rate of 
the fluid on the basis of the amplitude ratio obtained 
by the signal conversion section. 
Brief Description of the Drawings 

Fig. 1 is a view for explaining the basic 

25 principle of an electromagnetic flowmeter according to 
the present invention; 

Fig. 2 is a view showing an eddy current and 
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interelectrode electromotive force when the flow rate of 
a fluid to be measured is 0; 

Fig. 3 is a view showing an eddy current and 
interelectrode electromotive force when the flow rate of 
5 a fluid to be measured is not 0; 

Fig. 4 is a block diagram showing the 
arrangement of an electromagnetic flowmeter according to 
the first embodiment of the present invention; 

Fig. 5 is a graph showing the complex vector 
10 of the frequency component of the carrier of the 
interelectrode electromotive force in the first 
embodiment of the present invention; 

Fig. 6 is a graph showing the complex vector 
of the frequency component of the sideband of the 
15 interelectrode electromotive force in the first 
embodiment of the present invention; 

Fig. 7 is a graph showing the complex vector 
of the frequency component of the sideband of an 
interelectrode electromotive force in the third 
20 embodiment of the present invention; 

Fig. 8 is a view showing another example of 
the exciting coil arrangement in the electromagnetic 
flowmeter according to the present invention ; 

Fig. 9 is a sectional view showing an example 
25 of the electrode used in the electromagnetic flowmeter 
according to the present invention; 

Fig. 10 is a sectional view showing another 
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example of the electrode used in the electromagnetic 
flowmeter according to the present invention; and 

Fig. 11 is a block diagram showing the 
arrangement of a conventional electromagnetic flowmeter. 
5 Description of the Preferred Embodiments 
[Basic Principle] 

Before a description of the basic principle of 
the present invention, generally known basic 
mathematical knowledge will be described. A cosine wave 
10 Acos(O)t) and sine wave Bsin(cot), which have the same 

frequency and different amplitudes, are synthesized into 
the following cosine wave. A and B are amplitudes, and 
co is an angular frequency. 

2 2 1/2 

Acos(cot) +Bsin(0)t) - ( A + B ) cos(cot - £) 
15 for 8 = tan^B/A) ...(1) 

To analyze the synthesis of equation (1), it 
is convenient to map the cosine wave Acos(cot) and sine 
wave Bsin(tot) onto a complex coordinate plane while 
plotting the amplitude A of the cosine wave Acos(cat) 
20 along the real axis and the amplitude B of the sine wave 
Bsin(cot) along the imaginary axis. 

More specifically, on the complex coordinate 

2 2 1/2 

plane, a distance (A + B ) from the origin gives the 
amplitude of the synthetic wave, and an angle £ = 
25 tan^B/A) with respect to the real axis gives the phase 
difference between the synthetic wave and 0)t. 

In addition, on the complex coordinate plane. 
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the following relation holds 

Cexp = Ccos( a ) + jCsin( e ) ...(2) 

Equation (2) is an expression of a complex 
vector. In equation (2), j is the imaginary number 
5 unit, C is the length of the complex vector, and 6 is 
the direction of the complex vector. Hence, to analyze 
the geometrical relationship on the complex coordinate 
plane, it is convenient to use conversion to a complex 
vector. 

10 In the following description, to explain a 

behavior that is exhibited by an interelectrode 
electromotive force and the manner the present invention 
uses the behavior, mapping to the complex coordinate 
plane and geometrical analysis using a complex vector 

15 are employed. 

First, an interelectrode electromotive force 
which is irrelevant to the flow rate (flow velocity) of 
a fluid to be measured per unit time will be described. 
As shown in Fig. 1, an electromagnetic flowmeter 

20 comprises a measuring pipe 1 through which a fluid to be 
measured flows and a pair of electrodes 2a and 2b which 
oppose each other in the measuring pipe 1 to be 
perpendicular to both the magnetic field applied to the 
fluid to be measured and an axis PAX of the measuring 

25 pipe 1 and also come into contact with the fluid to be 
measured. The electrodes 2a and 2b detect an 
electromotive force generated by the magnetic field and 
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the flow of the fluid to be measured. 

The electromagnetic flowmeter also comprises a 
first exciting coil 3a and second exciting coil 3b. In 
the electromagnetic flowmeter, a plane PLN which 
5 includes the electrodes 2a and 2b and is perpendicular 
to the direction of the measuring pipe axis PAX is 
defined as a boundary in the measuring pipe 1. At this 
time, the first exciting coil 3a and second exciting 
coil 3b apply asymmetrical magnetic fields to the fluid 
10 to be measured on both sides of the plane PLN, i.e., the 
boundary in the measuring pipe 1. 

Of the magnetic field generated from the first 
exciting coil 3a, a magnetic field component (magnetic 
flux density) Bl which is perpendicular to both an 
15 electrode axis EAX that connects the electrodes 2a and 

2b and the measuring pipe axis PAX on the electrode axis 
EAX, and of the magnetic field generated from the second 
exciting coil 3b, a magnetic field component (magnetic 
flux density) B2 which is perpendicular to both the 
20 electrode axis EAX and the measuring pipe axis PAX on 
the electrode axis EAX are given by 

Bl = blcos(0)0t - 01) ....(3) 
B2 = b2cos(C00t - 92) ...(4) 
In equations (3) and (4), bl and b2 are the 
25 amplitudes, 0)0 is the angular frequency, and Ol and 0 2 
are the phase differences (phase delays) from coot. The 
magnetic flux density Bl will be referred to as the 
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magnetic field Bl, and the magnetic flux density B2 will 
be referred to as the magnetic field B2. 

An electromotive force caused by a change in 
magnetic field is obtained by a time differential dB/dt 
5 of the magnetic field. The magnetic fields Bl and B2 
generated from the first exciting coil 3a and second 
exciting coil 3b are differentiated as follows. 

dBl/dt = -blC0 0sin( 0)0t - 9 1) ...(5) 
dB2/dt = -b2C0 0sin( coot - 9 2) ...(6) 
10 When the flow rate of the fluid to be measured 

is 0, eddy currents generated by the magnetic fields Bl 
and B2 contain only components generated by a change in 
magnetic fields. An eddy current la by the magnetic 
field Bl and eddy current lb by the magnetic field B2 
15 have directions as shown in Fig. 2. 

Hence, in the plane that includes the 
electrode axis EAX and measuring pipe axis PAX, an 
interelectrode electromotive force Ea that is generated 
by a change in magnetic field Bl and is irrelevant to 
20 the flow rate (flow velocity) and an interelectrode 

electromotive force Eb that is generated by a change in 
magnetic field B2 and is irrelevant to the flow rate 
(flow velocity) have opposite directions, as shown in 
Fig. 2. 

25 At this time, a total interelectrode 

electromotive force E obtained by adding the 
interelectrode electromotive forces Ea and Be 
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corresponds to a value obtained by calculating the 
difference between the time differentials dBl/dt and 
dB2/dt of the magnetic fields and multiplying the 
difference by a coefficient k (a complex number related 
5 to the conductivity and dielectric constant of the fluid 
to be measured and the structure of the measuring pipe 
1). 

E = k{-b2C00sin( coot - 9 2) + blC00sin( COOt - 9l)} 

...(7) 

10 Equation (7) can be rewritten to 

E = -kb2C0 0sin( to Ot) cos ( - 9 2) 

- kb2C00cos( C00t)sin(- 9 2) 

+ kblC00sin( co0t)cos( - 6 1) 

+ kblco0cos( coot)sin( - 9 1) 
15 = {-b2sin(- 9 2) 

+ blsin( - 9 1) }0)0kcos( coot) 

+ {-b2cos(- 9 2) 

+ blcos(-£l)}co0ksin(co0t) 

.-.(8) 

20 When equation (8) is mapped onto a complex 

coordinate plane based on coot, a real axis component Ex 
and imaginary axis component Ey are given by 

Ex = {-b2sin(-02) + blsin ( - 9 1 ) } CO Ok ...(9) 
Ey = {-b2cos(- 9 2) + blcos ( - 9 1 ) } co Ok ...(10) 
25 Ex and Ey in equations (9) and (10) are 

rewritten to 

Ex = {-b2sin(-#2) + blsin( - 9 1 ) }O>0k 
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= {-b2cos(7T/2 + 02) + blcos(7T/2 + 01)}COOk 
= {b2cos(-7t/2 + 02) + blcos(7t/2 + 01)}O)Ok 

...(11) 

Ey = {-b2cos(-02) + blcos( - 9 1 ) } cook 
5 = £-b2sin(7t/2 + 02) + blsin(7t/2 + 01)}coOk 

= £b2sin(-7T/2 + 02) + blsin(7t/2 + 01)}COOk 

...(12) 

to obtain a complex vector Ec given by 
Ec = Ex + jEy 
10 = {b2cos(-7t/2 + 92) 

+ blcos(7t/2 + 01)}<OOk 
+ j{b2sin(-7T/2 + 02) 
+ blsin( 7C/2 + 0 l)}C00k 
= {blcos(7t/2 + 01) 
15 + ;jblsin(7t/2 + 01)}coOk 

+ {b2cos(-7t/2 + 02) 
+ jb2sin( - 7t/2 + 02)}coOk 
= blCO0kexp{3( 71/ 2 + 0 1)} 
+ b2CO0kexp{ j(-7t/2 + 0 2)} 
20 ...(13) 
The above -described coefficient k can be 
converted into a complex vector given by 
k = rkcos(0OO) + jrksin(0OO) 

= rkexp( j 0 00) ... (14) 

25 In equation (14), rk is a proportional 

coefficient, and 0 00 is the angle of the vector k with 
respect to the real axis. The angle 0 00 changes in 
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accordance with a delay of the magnetic field with 
respect to the exciting current or a change in 
conductivity of the fluid. The change in angle 0 00 is 
a flow rate measurement error. 

When equation (14) is substituted into 
equation (13), the interelectrode electromotive force Ec 
(an interelectrode electromotive force which is caused 
only by a time-rate change in magnetic field and is 
irrelevant to the flow velocity) converted into the 
complex vector is given by 

Ec = blcoOkexp( j( 7t/2 + 0 1)} 

+ b20)0kexp( j(-7C/2 + 02)} 
= blw0rkexp(j(7t/2 + 0 1 + 0 00)} 
+ b20)0rkexp(j(-7t/2 + 0 2 + 000)} 

. . . (15) 

in equation (15), bio) 0rkexp{ j ( K /2 +01+0 
00)} is a complex vector whose length is blcoOrk and 
angle from the real axis is */2 + 0 1 * 000, and b2C0 
0rkexp{ ;)(-*/ 2 + 0 2 + 000)} is a complex vector whose 
20 length is b2CO0rk and angle from the real axis is -7C/2 

+ 0 2 + 0 00. 

The interelectrode electromotive force caused 

by the flow rate (flow velocity) of the fluid to be 

me asured will be described next. When the flow velocity 

25 of the fluid to be measured is V (V * 0), eddy currents 

by the magne tic fields Bl and B2 respectively contain 

components V X Bl and V X B2 caused by the flow 
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velocity in addition to the eddy current components la 
and lb when the flow velocity is 0. For this reason, an 
eddy current la' by the magnetic field Bl and an eddy 
current lb 1 by the magnetic field B2 have directions as 
5 shown in Fig. 3. 

Hence, an interelectrode electromotive force 
Ea' generated by the flow velocity V of the fluid to be 
measured and the magnetic field Bl and an interelectrode 
electromotive force Eb 1 generated by the flow velocity V 

10 and the magnetic field B2 have the same direction. 

At this time, a total interelectrode 
electromotive force Ev obtain by adding the 
interelectrode electromotive forces Ea' and Eb 1 
generated by the flow velocity corresponds to the sum of 

15 a value obtained by multiplying the magnetic field Bl by 
a coefficient kv (a complex number related to the flow 
velocity V, the conductivity and dielectric constant of 
the fluid to be measured, and the structure of the 
measuring pipe 1) and a value obtained by multiplying 

20 the magnetic field B2 by the coefficient kv. 

Ev - kv{blcos(C00t - Ol) + b2cos(C00t - #2)} 

... (16) 

When the term of sin and the term of cos of 
equation (16) are expanded, we obtain 
25 Ev = kvblcos( wot ) cos (- 9 1) 

- kvblsin(co0t)sin(-(9l) 
+ kvb2cos( coot )cos( - 6 2) 
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- kvb2sin(C00t)sin(-#2) 
= {blcos(- 01) 

+ b2cos( - 9 2) }kvcos( coot) 
+ {-blsin(- 9 1) 

5 - b2sin(-0 2)}kvsin(coot) ...(17) 

When equation (17) is mapped onto the complex 
coordinate plane based on coot, a real axis component 
Evx and imaginary axis component Evy are given by 

Evx = {blcos(-01) + b2cos(- 9 2)}kv ...(18) 
10 Evy = {-blsin(- (9 1) - b2sin( - 9 2 ) }kv ...(19) 

Equations (18) and (19) are transformed into a 
complex vector Eve. 

Evx = {blcos(-fll) + b2cos(-#2)}kv 

= {blcos(#l) + b2cos( 0 2)}kv ...(20) 
15 Evy = {-blsin(-#l) - b2sin( - 9 2 ) }kv 

= {blsin(<9l) + b2sin(0 2)}kv ...(21) 

Eve = Evx + jEvy 

= {blcos(01) + b2cos(0 2)}kv 
20 + j{blsin(01) + b2sin(0 2)}kv 

= {blcos(/9l) + jblsin(#l)}kv 

+ {b2cos(0 2) + jb2sin( 9 2)}kv 
= blkvexp(j 0 1) + b2kvexp(j (9 2) 
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25 The above-described coefficient kv is 

transformed to a complex vector. 

kv = rkvcos(<9 01) + jrkvsin( 9 01 ) } 
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= rkvexp( 3 9 01) ... (23) 

In equation (23), rkv is a proportional 
coefficient, 0 01 is the angle of the vector kv with 
respect to the real axis. In this case, rkv corresponds 
5 to a value obtained by multiplying the proportional 
coefficient rk (equation (14)) by the flow velocity V 
and a proportional coefficient 7, so v = VT . That is, 
rkv = rkVT ... (24) 

When equation (23) is substituted into 
10 equation (22), the interelectrode electromotive force 
Eve converted into complex coordinates is obtained as 
Eve = blkvexp(j01) + b2kvexp(j0 2) 
= blrkvexp{ ( j ( 0 1 + 0 01)} 

+ b2rkvexp{ (j( 0 2 +001)} ...(25) 
15 In equation (25), blrkvexp{ j ( 9 1 + 001)} is a 

complex vector whose length is blrkv and angle from the 
real axis is 01 + 0 01, and b2rkvexp{ j ( 0 2 + 0 01)} is a 
complex vector whose length is b2rkv and angle from the 
real axis is 0 2 + 0 01. 
20 From equations (15) and (25), a total 

interelectrode electromotive force Eac obtained by 
adding the interelectrode electromotive force Ec 
generated by a time-rate change in magnetic field and 
the interelectrode electromotive force Eve generated by 
25 the flow velocity of the fluid is given by 
Eac - Ec + Eve 

= blC0 0rkexp{ j ( 7Z/2 + 01 + 0 00)} 



- 16 - 



+ b2 CO Orkexp{ j ( - 71/2 + 0 2 + 0 00)} 
+ blrkvexp{ j( 0 1 + 0 01)} 

+ b2rkvexp{ j ( 0 2 +001)} ...(26) 
As is apparent from equation (26), the 
5 interelectrode electromotive force Eac is described by 
the four complex vectors blo>0rkexp{ j( re 1 2 +01+0 

00) }, b20)0rkexp{j( - 7T/2 + 0 2 + 0 00)}, blrkvexp{ j ( 0 1 + 

001) }, and b2rkvexp{ j ( 0 2 + 001)}. 

The length of a synthetic vector obtained by 

10 synthesizing the four complex vectors represents the 
amplitude of the output (interelectrode electromotive 
force Eac), and an angle <S> of the synthetic vector 
represents the phase difference (phase delay) of the 
interelectrode electromotive force Eac from the phase co 

15 Ot of the input (exciting current). 

In the present invention, the carrier of an 
angular frequency 0)0 is amplitude-, phase-, or 
frequency-modulated by the modulated wave of an angular 
frequency (0 2 to obtain an exciting current. The 

20 exciting current is supplied to the first and second 

exciting coils 3a and 3b to apply asymmetrical magnetic 
fields to the fluid to be measured on both sides of the 
plane PLN, i.e., the boundary in the measuring pipe 1. 
Accordingly, a plurality of frequency 

25 components coo, COO + £o>2, and co 0 - £ 0)2 (£ is an 

integer (£ ^ 1); for amplitude modulation, only I = 1) 
are generated in the interelectrode electromotive force 
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Eac. From two of these frequency components, an 
asymmetrical exciting characteristic parameter 
(amplitude ratio or phase difference) is obtained, which 
depends on the flow rate of the fluid and does not 
5 depend on the variation in delay (000) of the magnetic 
field with respect to the exciting current or the shift 
of the amplitude of the magnetic field. On the basis of 
the asymmetrical exciting characteristic parameter, a 
flow rate measurement error due to the variation in 

10 delay of the magnetic field with respect to the exciting 
current or the shift of the amplitude of the magnetic 
field is automatically corrected. This is the basic 
technical idea of the present invention. 

Accordingly, in-phase component noise is 

15 removed so that the rectangular wave exciting method 

need not be used, and the sine wave exciting method can 
be used. 

[First Embodiment] 

An embodiment of the present invention will 

20 be described below in detail. Fig. 4 shows the 

arrangement of an electromagnetic flowmeter according to 
the first embodiment of the present invention. The same 
reference numerals as in Fig. 1 denote the same 
components in Fig. 4. 

25 The electromagnetic flowmeter according to 

this embodiment comprises a measuring pipe 1 , electrodes 
2a and 2b, first and second exciting coils 3a and 3b, 
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and a power supply section 4 which supplies a first 
exciting current to the first exciting coil 3a and a 
second exciting current to the second exciting coil 3b. 

The electromagnetic flowmeter also comprises a 
5 signal conversion section 5 and flow rate output section 
6 . The signal conversion section 5 obtains an amplitude 
by separating a component having an angular frequency co 
0 from an electromotive force detected by the electrodes 
2a and 2b, obtains an amplitude by separating a 

10 component having a sum frequency (coo + (0 2) or a 
difference frequency (coO - 00 2) of the angular 
frequency coo and an angular frequency CO 2 from the 
electromotive force, and obtains the ratio of these 
amplitudes . The flow rate output section 6 calculates 

15 the flow rate of a fluid to be measured on the basis of 
the amplitude ratio obtained by the signal conversion 
section 5, 

The first exciting coil 3a is arranged 
downstream of a plane PLN at a position separated from 

20 it by an offset distance dl . The second exciting coil 
3b is arranged upstream of the plane PLN at a position 
separated from it by an offset distance d2 . That is, 
the second exciting coil 3b is arranged on the opposite 
side of the first exciting coil 3a with respect to the 

25 plane PLN 

The power supply section 4 supplies a first 
sine wave exciting current having the first angular 
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frequency coo to the first exciting coil 3a. In this 
embodiment, bl = b and 9 1 = 0 in equation (3). Of the 
magnetic field generated from the first exciting coil 3a 
when the first exciting current is supplied from the 
5 power supply section 4, a magnetic field component Bl 
that is perpendicular to both an electrode axis EAX and 
a measuring pipe axis PAX on the electrode axis EAX is 
given by 

Bl = bcos( coot) ... (27) 

10 The power supply section 4 also supplies a 

second exciting current to the second exciting coil 3b. 
The second exciting current is obtained by 
amplitude-modulating a sine wave carrier having the same 
angular frequency coo as that of the carrier component 
15 of the first exciting current and a predetermined phase 
difference 6 2 by using a modulated sine wave having the 
second angular frequency co 2 . 

Of the magnetic field generated from the 
second exciting coil 3b when the second exciting current 
20 is supplied from the power supply section 4, an 

amplitude b2 of a magnetic field component B2 that is 
perpendicular to both the electrode axis EAX and the 
measuring pipe axis PAX on the electrode axis EAX is 
given by 

25 b2 = b{l + m a cos(C02t)} ...(28) 

In equation (28), m„ is an amplitude 
modulation index. From equations (4) and (28), the 
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magnetic field B2 is given by 

B2 = Ml * m cos (0>2t)}cos( coot - 02) ■-■^9) 
in equation (26). bl - b. 01 = 0, and 001 = 
0 00. When the magnetic fields Bl and B2 are given by 
equations (27) and (29), we obtain 
Eac = Ec + Eve 

= bCO0rkexp{3(7t/2 + 0 00)} 

+ b{ i + mcos(C02t)}C00rkexp{j(-7C/2 
+ 9 2 + 0 00)} + brkvexp{j( 000)} 
+ b {l + mcos(co2t)}rkvexp{j(0 2 + 000)} 
= bco0rkexp{j(7t/2 + 000)} 

+ bC00rkexp{j(-7t/2 + 0 2 + 000)} 

+ brkvexp{ j( 0 00)} 

+ brkvexp{j( 0 2 + 0 00)} 

+ mC os(C02t)bC00rkexp{3(-^/2 + ^ + 000)} 

15 a 

+ m cos(w2t)brkvexpU< 02 * 600)} 

...(30) 

Four vectors on the right-hand side of 
equation (30). i.e.. bO>0rKe*p( ) ( */» * "00,) - the 
20 first term. bO>0r*e*p< 3 ( - */2 * "2 ♦ * s the 

second tern,. brkvexpU ( « 00 ) } as the third tern,, and 
orKvex P U<92 ♦ «00» as the fourth tern, correspond to 
fundamental vectors obtained when no amplitude 

modulation is done. 

The vector as the fifth term on the right-hand 
„ nnl i e mcos(C02t)b(00rkexp{j(-7C/2 

S1UC . 

+ 9 2 + 0 00)} can be rewritten to bCOOrkm cost coot - (-* 
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/2 + 0 2 + 0 00) }cos( C0 2t ) as time expression. This time 
expression can further be rewritten to 
bcoOrkm cos { coot 

a 

- (-7T/2 + 02 + 0OO)}cos(CO2t) 
5 = ( 1/2 )bco0rkm a [cos{ coot 

- (-7T/2 + 02 + 000) + C02t} 

+ cos{C00t - (-7C/2 + 02 + 000) - C02t}] 
= (l/2)bC00rkm a cos{(CO0 + C02)t 

- (-7T/2 + 0 2 + (9 00 ) } 

10 + (l/2)bC00rkm a cos{ ( COO + C02)t 

- (-7T/2 + 02 + 000)} 

... (31) 

As is apparent from equation (31), the fifth 
term on the right-hand side of equation (30) forms a 
15 vector ( 1/2 )bco 0rkm.exp{ j ( - % /2 + 0 2 + 0 00)} on each of 
a complex plane based on the angular frequency (coo + co 
2) and a complex plane based on the angular frequency 
(coo - C02) . 

The vector as the sixth term on the right-hand 
20 side of equation (30), i.e., mcos( C02t )brkvexp{ j ( 0 2 + 
000)} can be rewritten to brkvm a cos{ coot - (02 + 0 
00) }cos( co2t ) as time expression. This time expression 
can further be rewritten to 

brkvm a cos{ COOt - ( 02 + 0 00 ) }cos ( CO 2t ) 
25 = ( 1/2 )brkvm a tcos{ coot - (02 + 0 00) + C02t} 

+ cos{C00t - ( 0 2 + 000) - C02t}] 
= (l/2)brkvm a cos{(C0 0 + C02)t - (02 + 0 00)} 
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+ l/2brkvm a cos{ ( coo - co2)t 
- ( 02 + 000)} 

. . . (32) 

As is apparent from equation (32), the sixth 
5 term on the right-hand side of equation (30) forms a 
vector ( l/2)brkvm a exp{ j ( 9 2 + 0 00)} on each of the 
complex plane based on the angular frequency (0)0 + 0)2) 
and the complex plane based on the angular frequency ( to 
0 - CO 2) . 

10 As is apparent from the above description, the 

fifth and sixth terms on the right-hand side of equation 
(30) form a complex vector Earn on each of the complex 
plane based on the angular frequency (COO + 0J2) and the 
complex plane based on the angular frequency (coo - CO 
15 2). The complex vector Earn is given by 

Earn = (l/2)bco0rkm a exp{ j(-7t/2 + 02 + 000)} 
+ (l/2)brkvm a exp{j( 0 2 + 0 00)} 

. . • (33) 

In addition, the first to fourth terms on the 
20 right-hand side of equation (30) form a complex vector 
Eor on a complex plane based on the angular frequency co 
0. The complex vector Eor is given by 
Eor = bC0 0rkexp{j( 7T/2 + 0 00)} 

+ bCO0rkexp{j(-7C/2 + 02 + 0 00)} 
25 + brkvexp{ j ( 9 00) } 

+ brkvexp{ j ( 0 2 + 0 00 ) } 

. . . (34) 
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Fig. 5 shows the complex vector Eor formed on 
the complex plane based on the angular frequency coo of 
the carrier. Fig. 6 shows the complex vector Earn formed 
on the complex plane based on the angular frequency (0)0 
5 + C02)or(co0-co2)of the sideband. 

The complex vector Earn is the synthetic vector 
of two complex vectors (l/2)b0)0rkinexp{ j( - K/2 +02+0 
00)} and ( 1/2 )brkvinexp{ j ( 0 2 + 0 00)}, which are 
perpendicular to each other. This synthetic vector only 

10 rotates about an origin A of the complex coordinate 
system shown in Fig. 6 as the angle 9 00 changes. 
Hence , the synthetic vector has a predetermined 
magnitude independently of the change in angle 0 00. 
The magnitude |Eam| of the complex vector Earn is given 

15 by 

| Earn | - [{(l/2)bcoOrkm a } 2 

+ {(l/2)brkvm a } 2 ] 1/2 

2 2 1/2 

= (l/2)brkm a {0)0 + v ) 

...(35) 

20 The complex vector Eor is a synthetic vector 

of four complex vectors ba>0rkexp{ j( it 12 + 0 00)}, ba> 
0rkexp{ j(-7T/2 + 0 2 + 000)}, brkvexp{ j ( 0 00 ) } , and 
brkvexp{j ( 0 2 + 000)}, which form two isosceles 
triangles AABC and ACDE, as shown in Fig. 5. This 

25 synthetic vector only rotates about an origin B of the 
complex coordinate system shown as the angle 0 00 
changes. Hence, the synthetic vector has a 
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predetermined magnitude independently of the change in 
angle #00. When the bases of the isosceles triangles 
AABC and AcDE are synthesized, the magnitude |Eor| of 
the complex vector Eor is given by 
5 I Eor | = 2b0)0rksin( 9 2/2) 

+ 2brkvcos( 0 2/2) 
= 2brk{o>0sin( 0 2/2) + vcos(02/2)} 

... (36) 

A ratio Ram of the magnitude |Eor| of the 
10 complex vector Eor to the magnitude |Eam| of the complex 
vector Earn is given by 
Ram = | Eor | / | Earn) 

= [2brk{C00sin( 9 2/2) + vcos(02/2)}] 
/{(l/2)brkm a {<W + vY /2 } 
15 = [4{0)0sin( 92/2) + vcos(02/2)}] 

2 2 1/2 

/{m a {co0 + v ) } 

. . . (37) 

The angular frequency too, phase difference 0 
2, and amplitude modulation index m a are irrelevant to 
20 the angle 0 00 or the amplitude b of the magnetic field 
Bl (the amplitude of the carrier component of the 
magnetic field B2). Equation (37) that represents the 
ratio Ram has no term containing the angle 9 00 or 
amplitude b. 

25 Hence, even when the angle 0 00 changes, or 

the amplitude b shifts, the ratio Ram does not change. 
When the magnitude |Eor| of the complex vector Eor and 
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the magnitude |Eam| of the complex vector Earn are 
obtained, and the flow rate is detected on the basis of 
the ratio Ram, a flow rate measurement error due to the 
variation in delay of the magnetic field with respect to 
5 the exciting current or the shift of the amplitude of 
the magnetic field can be quickly and automatically 
canceled. 

To obtain the flow rate of the fluid to be 
measured, equation (37) is rewritten to 

2 2 1/2 

10 v = 0)0{-8sin( 6 2) + Ramm a (16 - RammJ } 

2 2 

/{8 + 8cos(0 2) - Ramm a } 

..(38) 

From equation (24), equation (38) can be 
rewritten to 

si 2 2 1/2 

15 V = a X 0)0{-8sin( 9 2) + Ramm a (16 - RammJ } 

2 2 

/{8 + 8cos(0 2) - Ramm } 
for a = 1/ T ... (39) 

where a (or 7 ) is a predetermined coefficient. The 
signal conversion section 5 detects the electromotive 

20 force Eac between the electrodes 2a and 2b and 
frequency- separates the detected interelectrode 
electromotive force Eac by a filter. Accordingly, the 
amplitude of the component having the angular frequency 
(0)0 + 0)2) or (COO - 0)2) (the magnitude |Eam| of the 

25 complex vector Earn) and the amplitude of the component 
having the angular frequency 0)0 (the magnitude |Eor| of 
the complex vector Eor) are obtained. 
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The signal conversion section 5 calculates the 
ratio Ram of the magnitudes |Eor| to |Eam|. 

The interelectrode electromotive force Eac can 
also be frequency- separated even by using a bandpass 
5 filter. Actually, when a comb- shaped digital filter 
called a comb filter is used, the interelectrode 
electromotive force can easily be separated into the 
three frequency components coo, (coo + C02), and ( <x> 
0 - 0)2) . 

10 The flow rate output section 6 calculates a 

flow velocity V of the fluid to be measured, i.e., the 
flow rate per unit time by using equation (39) on the 
basis of the ratio Ram obtained by the signal conversion 
section 5. With the above arrangement, in this 

15 embodiment, the flow rate can be calculated while 
quickly and automatically canceling the flow rate 
measurement error so that accurate flow rate measurement 
can be executed. 
[Second Embodiment] 

20 The second embodiment of the present invention 

will be described next . The arrangement of an 
electromagnetic flowmeter according to this embodiment 
is the same as in the first embodiment and will be 
described with reference to Fig. 4. The operation of a 

25 power supply section 4 according to this embodiment is 
the same as in the first embodiment. A magnetic field 
Bl generated from a first exciting coil 3a when an 
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exciting current is supplied from the power supply 
section 4 is given by equation (27). A magnetic field 
B2 generated from a second exciting coil 3b is given by 
equation (29). When bl = b, 01 = 0, and 9 01 = 9 00 in 
5 equation (26), an interelectrode electromotive force Eac 
is given by equation (30). 

Four vectors that represent a component 
(complex vector Eor) having an angular frequency WO of 
the interelectrode electromotive force Eac have a 

10 geometrical relationship as shown in Fig. 5 described in 
the first embodiment. Two vectors that represent a 
component (complex vector Earn) having an angular 
frequency (coo + 0)2) or ( co 0 - C02) have a geometrical 
relationship as shown in Fig. 6. 

15 From Fig. 5, a phase difference <f>or between 

the complex vector Eor and the first exciting current 
supplied to the first exciting coil 3a is given by 

</>or = 9 2/2 + 9 00 ... (40) 

From equation (40), an angle 9 00 is given by 

20 9 00 = <f>or - 9 2/2 . . . (41) 

Let <£am be the phase difference between the 
component having the angular frequency (coo + 0)2) of 
the interelectrode electromotive force Eac and the 
component having the angular frequency (0)0 + (0 2) when 

2 5 02 = 0 in the second exciting current supplied to the 
second exciting coil 3b or the phase difference between 
the component having the angular frequency (0)0 - 0)2) 
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of the interelectrode electromotive force Eac and the 
component having the angular frequency ( co 0 - <j) 2 ) when 
9 2 = 0 in the second exciting current. 

In addition, let <t>x be the angle made by a 
5 vector ( 1/2 )b0) Orkm & exp{ j ( - % 12 + 0 2 + 0 00)} as the 

first term on the right-hand side of equation (33) and a 
vector ( l/2)brkvm a exp{ j( 9 2 + 000)} as the second term 
on the right-hand side. 

At this time, a relation given by equation 
10 (42) holds between the phase difference <£>am and the 
angle <t>x as is apparent from Fig. 6. In the example 
shown in Fig. 6, however, the phase difference <f>am has 
a negative value. 

0x = 7T/2+ <£am - 0 00- 02 
15 = %/2 + (£am - (<i>or - 02/2) - 02 

= 71/2 + <£>am - <f>or - 02/2 ...(42) 
In addition, AABC shown in Fig. 6, which is 
formed by the two complex vectors ( 1/2 )bco 0rkinexp{ j ( - K 
/2 + 0 2 + 0 00)} and ( 1/2 )brkvinexp{ j ( 9 2 + 000)} and 
20 their synthetic vector Earn, is a right triangle. Hence, 
tan<£x = {(l/2)brkvm a }/{(l/2)bco0r3cm a } 

= v/ 000 ... (43) 

Equation (43) can be rewritten to 
v = 0)0tan(f>x = Otf0tan{ 7T/2 + (/> am - </)or - 02/2} 
25 ...(44) 
From equation (24), equation (44) can be 
rewritten to 
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V = a X C00tan(7T/2 + 4>am - <bor - 92/2) 
for a = 1/ r ... (45) 

Equation (45) that represents a flow velocity 
V has no terms of the amplitude of the magnetic field Bl 
5 (the amplitude of the carrier component of the magnetic 
field B2) and angle 0 00. Hence, when the flow rate is 
detected on the basis of the phase differences <l>am and 
<f>or, a flow rate measurement error due to the variation 
in delay of the magnetic field with respect to the 
10 exciting current or the shift of the amplitude of the 
magnetic field can be quickly and automatically 
canceled. 

A signal conversion section 5 
frequency- separates the interelectrode electromotive 

15 force Eac by using a filter, as in the first embodiment, 
to obtain the phase difference (f>or between the complex 
vector Eor and the first exciting current supplied to 
the first exciting coil 3a. The signal conversion 
section 5 also frequency- separates the second exciting 

20 current supplied to the second exciting coil 3b by using 
the filter to obtain the phase difference <£am between 
the complex vector Earn and the component having the 
angular frequency (coo + C02) or (coo - 0)2) of the 
second exciting current. 

25 A flow rate output section 6 calculates the 

flow velocity V of the fluid to be measured by using 
equation (45) on the basis of the phase differences $am 
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and <f>or obtained by the signal conversion section 5. 
With the above arrangement, the same effect as in the 
first embodiment can be obtained. 

Instead of using the component having the 
5 angular frequency (coo + CO 2) or (coo - 0)2) of the 

exciting current supplied to the second exciting coil 3b 
as the reference phase of the phase difference 0am, the 
phase difference 0 am may be detected with reference to 
a timing when cos(C0 2t) becomes 1.0. 

10 [Third Embodiment] 

The third embodiment of the present invention 
will be described next. The arrangement of an 
electromagnetic flowmeter according to this embodiment 
is the same as in the first embodiment and will be 

15 described with reference to Fig, 4. A power supply 

section 4 according to this embodiment supplies a first 
exciting current to a first exciting coil 3a. The first 
exciting current is obtained by amplitude -modulating a 
sine wave carrier having a first angular frequency 0)0 

20 by a modulated sine wave having a second angular 

frequency C0 2. Of the magnetic field generated from the 
first exciting coil 3a when the first exciting current 
is supplied from the power supply section 4, a magnetic 
field component Bl which is perpendicular to both an 

25 electrode axis EAX and a measuring pipe axis PAX on the 
electrode axis EAX has an amplitude bl given by 

bl = b{l - m a cos(C0 2t)} ...(46) 
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In this embodiment, 01 = 0 in equation (3). 
From equations (3) and (46), the magnetic field 
component Bl is given by 

Bl = b{l - m a cos( C0 2t)}cos( coot) ...(47) 
5 The power supply section 4 also supplies a 

second exciting current to a second exciting coil 3b. 
The second exciting current is obtained by 
amplitude-modulating a sine wave carrier having the same 
angular frequency coo as that of the carrier component 
10 of the first exciting current and a predetermined phase 
difference 9 2 by a modulated sine wave having the same 
angular frequency 0)2 as that of the modulated wave 
component of the first exciting current and an opposite 
phase. 

15 Of the magnetic field generated from the 

second exciting coil 3b when the second exciting current 
is supplied from the power supply section 4, a magnetic 
field component B2 which is perpendicular to both the 
electrode axis EAX and the measuring pipe axis PAX on 
20 the electrode axis EAX has an amplitude b2 given by 

b2 = b{l + m a cos(C0 2t)} ...(48) 
From equations (4) and (48), the magnetic 
field component B2 is given by 

B2 = b{l + m a cos( C0 2t) }cos( OOOt - 9 2) ...(49) 
25 In equation (26), 01 = 0, and 001 = 000. 

when the magnetic fields Bl and B2 are given by 
equations (47) and (49), we obtain 
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Eac = Ec + Eve 

= b {l - mcos(0)2t)}C0 0rkexp{j(7t/2 + ^00)} 
+ b{l + mcos(w2t)}C00rkexp{j(-7t/2 + 02 

+ 000)} 

+ b{l - mcos(co2t)}rkvexp{j(0OO)} 
+ b{l + mcos(co2t)}rkvexp{j(0 2 + 000)} 
= bC0 0rkexp{j(7C/2 + 000)} 

+ b 0)0rkexp{3(-^/2 + 0 2 + 0 00)} 

+ brkvexp{j( # 00)} 

+ brkvexp{j(0 2 + 000)} 

- mcos(w2t)b(00rkexp{j(7r/2 + #00)} 

- mcos( W2t)brkvexp{d( #00)} 

+ m a cos(co2t)bcoOrkexp{j(-7r/2 + 02 + 000)} 
+ m cos(O)2t)brkvexp{j(0 2 + 000)} 
= bWOrkexp{j(7t/2 + 0 00)} 

+ ba)0rkexp{j(-7t/2 + 02 + 000)} 

+ brkvexp{j( 0 00)} 

+ brkvexp{j( 0 2 + 000)} 

+ mcos((02t)bco0rkexp{j(-^/2 + ^00)} 

+ m cos(C0 2t)brkvexp{j(7t + 0 00)} 

+ m a cos(0>2t)bC00rkexp{D(-7r/2 + 0 2 + 000)} 

+ mcos(co2t)brkvexp{j(0 2 + 0 00)} 

...(50) 

Four vectors on the right-hand side of 
equation (50), i.e., bco 0rkexp( 3 ( K ,2 * 9 00)} as the 
£irst term< b C0 0rkexp{j(-7t/2 + 02 + 000)} as the 
secld term. brkvexpt 3 ( 0 00 ) } as the third term, and 
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b rKvesp«(o2 + 0OO)> as the fourth term correspond to 
fundamental vectors obtained when no amplitude 

modulation is done. 

The vector as the fifth term on the right -hand 
side of equation (50), i.e., mcos( co2t)bC0 0rkexp( j< - * /2 
+ 000)} can be rewritten to bcoorkm cos { coot - (-*/2 + 
0OO)}cos(co2t) as time expression. This time 
expression can further be rewritten to 

r ,. nf (-7Z/2 + 9 00) }cos( C02t) 
bco0rkm a cos{ coot - (-n / z + v w i j 

= (l/2)bC00rkm a [cos{C00t - (-*/2 + 9o °) 

+ C02t} + cos{C00t - (-7C/2 + 000) - C02t>] 
= (l/2)bC00rkm a cos{ ( C00 + co2)t 

- {-7Z/2 + 000)} 

+ (i/2)bC0 0rkm a cos{ ( coo + co2)t 

_ (-7t/2 + 000)} 

. . .(51) 

As is apparent from equation (51), the fifth 
term on the right-hand side of equation (50) forms a 
vector (l/2)bC00rkmexp{ j (- 7 r/2 * 000)} on each of a 
, complex plane based on an angular frequency (COO + co2) 
and a complex plane based on an angular frequency ( co 
0 - C02) . 

The vector as the sixth term on the right-hand 

Knl i e mcos(C0 2t)brkvexp{ j( 7C + 6 
side of equation (50), i.e., iu^=>\ 

5 00)} can be rewritten to brkvm cos( co 0t - ( it + B 

„, time expression. This time expression 

can further be rewritten to 
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brkvm a cos{Q)Ot - ( it + 0 00 ) }cos ( co 2t ) 

= (l/2)brkvm [cos{C00t - ( it + 000) + 0)2t> 

+ cos{C00t - (7t + 000) - C02t}] 
= (l/2)brkvincos{ ( coo + co2)t - ( 7t + 000)} 
5 + (l/2)brkvmcos{(C00 + C02)t 

- ( 7t + 000)} 

... (52) 

As is apparent from equation (52), the sixth 
term on the right-hand side of equation (50) forms a 
10 vector ( 1/2 )brkvm a exp{ j ( it + 0 00)} on each of the 

complex plane based on the angular frequency (coo + 0)2) 
and the complex plane based on the angular frequency ( co 
0 - 0)2) . 

The seventh term on the right-hand side of 
15 equation (50), i.e., incos ( Q) 2t )btc>0rkexp{ J ( - % /2 + 02 + 
0 00)} is the same as the fifth term on the right-hand 
side of equation (30). Hence, the seventh term on the 
right-hand side of equation (50) forms a vector (l/2)bco 
0rkm g exp{ j ( - it /2 + 02 + 0 00)} on each of the complex 
20 plane based on the angular frequency (coo + co2) and the 
complex plane based on the angular frequency (coo - co 
2). 

The eighth term on the right-hand side of 
equation (50), i.e., mcos( co2t )brkvexp{ j ( 0 2 + 000)} is 
25 the same as the sixth term on the right-hand side of 

equation (30). Hence., the eighth term on the right-hand 
side of equation (50) forms a vector ( 1/2 )brkvinexp{ j ( 6 2 
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+ 0 00)} on each of the complex plane based on the 
angular frequency (coo + C0 2) and the complex plane 
based on the angular frequency (coo - 0)2), 

As is apparent from the above description, the 
5 fifth, sixth, seventh, and eighth terms on the 

right-hand side of equation (50) form a complex vector 
Earn on each of the complex plane based on the angular 
frequency (coo + (0 2) and the complex plane based on the 
angular frequency (0)0 - (0 2). The complex vector Earn 
10 is given by 

Earn = (l/2)bC00rkm a exp{ j(-7T/2 + 0 00)} 

+ (l/2)bCO0rkm a exp{ j(-7C/2 + 02 + 000)} 
+ (l/2)brkvm a exp{ j ( K + 0 00)} 
+ (l/2)brkvm a exp{ j( 0 2 + 0 00)} 
15 ...(53) 
In addition, the first to fourth terms on the 
right-hand side of equation (50) form a complex vector 
Eor given by equation (34) on a complex plane based on 
the angular frequency coo. Four vectors that represent 
20 the complex vector Eor have a geometrical relationship 
as shown in Fig. 5 described in the first embodiment. 

Four vectors that represent the complex vector 
Earn, i.e., ( 1/2 )bco 0rkm a exp{ j ( - 7t 12 + 0 00)}, (l/2)b0) 
0rkm a exp{j(-7r/2 + 0 2 + 0 00)}, ( 1/2 )brkvinexp{ j ( K + 0 
25 00)}, and ( 1/2 )brkvm.exp{ j ( 0 2 + 0 00)} have a 
geometrical relationship as shown in Fig. 7. 

As shown in Fig. 7, the complex vector Earn is 
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the synthetic vector of four vector (l/2)bco 
Orkm a exp{ j(-7T/2 + 0 00)}, ( 1/2 )b0) 0rkmexp{ j ( - K /2 + 9 2 + 
0 00)}, (l/2)brkvmexp{ j ( K + 9 00)}, and 
( 1/2 )brkvm a exp{ j ( 9 2 + 000)}, which form two isosceles 
5 triangles AABC and ACDE as shown in Fig. 7. This 

synthetic vector only rotates about an origin B of the 
complex coordinate system as the angle 0 00 changes. 
Hence, the synthetic vector has a predetermined 
magnitude independently of the change in angle 9 00 . 

10 When the bases of the isosceles triangles AABC and A 
CDE are synthesized, the magnitude |Eam| of the complex 
vector Earn is given by 

| Earn | = 2(l/2)bC00rkm a cos( 9 2/2) 
- 2(l/2)brkvm a sin( 9 2/2) 

15 = brkm a {C00cos( 9 2/2) + vsin(0 2/2)} 

. . . (54) 

The magnitude |Eor| of the complex vector Eor 
is given by equation (36). A ratio Ram of the magnitude 
| Eor | of the complex vector Eor to the magnitude |Eam| 
20 of the complex vector Earn, which is given by equation 
( 54 ) , is given by 

Ram = | Eor | / | Eam| 

= [2brk{W0sin( 9 2/2) + vcos(0 2/2)}] 
/[brkm a {0)0cos( 9 2/2) - vsin( (9 2/2 ) } ] 
25 = [2{CO0sin( 9 2/2) + vcos ( 9 2/2 ) } ] 

/[m a {C00cos( 9 2/2) - vsin( <9 2/2) }] 

. . . (55) 
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As in the first embodiment, equation (55) that 
represents the ratio Ram has no term containing the 
angle 9 00 or amplitude b. Hence, even when the angle 
9 00 changes, or the amplitude b shifts, the ratio Ram 
5 does not change. When the flow rate is detected on the 
basis of the ratio Ram, a flow rate measurement error 
due to the variation in delay of the magnetic field with 
respect to the exciting current or the shift of the 
amplitude of the magnetic field can be quickly and 
10 automatically canceled. 

To obtain the flow rate of the fluid to be 
measured, equation (55) is rewritten to 

v = 0)0{Ramm a cos( 9 2/2) - 2sin(<9 2/2)} 

/{Raram a sin( 9 2/2) + 2cos(#2/2)} ..(56) 
15 From equation (24), equation (56) can be 

rewritten to 

V= a X C00(Ramm a cos( 9 2/2) - 2sin(#2/2)} 

/{Ramm a sin( 9 2/2) + 2cos(0 2/2)} 
for a = 1/ r , . . (57) 

20 A signal conversion section 5 

frequency- separates the interelectrode electromotive 
force Eac by using a filter, as in the first embodiment, 
to obtain the amplitude of the component having the 
angular frequency (0)0 + (0 2) or ((0 0 - (0 2) (the 
25 magnitude |Eam| of the complex vector Earn). The signal 
conversion section 5 also obtains the amplitude of the 
component having the angular frequency (0 0 (the 
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magnitude |Eor| of the complex vector Eor, . The signal 
conversion section 5 calculates the ratio Ram of the 

magnitudes |Bor| to |Eam|. 

X flow rate output section 6 calculates a flow 
velocity V of the fluid to he measured h y usin g elation 
,57, on the hasis of the ratio Ram obtained hy the 
signal conversion section 5. With the above 

as in the first embodiment 
arrangement, the same effect as in 

can be obtained. 
[Fourth Embodiment] 

Th e fourth embodiment of the present invention 

will be described next. The arrangement of an 

electromagnetic flowmeter according to this embodiment 

=,«= in the first embodiment and will be 
is the same as in tne 

5 described with reference to Fig. 4. A power supply 

action 4 according to this embodiment supplies a first 
sine wave exciting current having a first angular 

,>n 4-0 a first exciting coil 3a. In this 
frequency 0)0 to a tirst 

6 mbodiment. bl - b. and 91 = 0 in eguation (3,. Of the 

2 0 magnetic field generated from the first exciting coll 3a 

when the first exciting current is supplied from the 

power supply section 4. a magnetic field component .1 

wmch is perpendicular to both an electrode axis RAX and 

. np axis pax on the electrode axis EAX is 
a measuring pipe axis w 

25 given by equation (27). 

The power supply section 4 also supplies a 
second exciting current to a second exciting coil 3b. 
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The second exciting current is obtained by 
phase -modulating a sine wave carrier having the same 
angular frequency 0)0 as that of the first exciting 
current and a predetermined phase difference 0 2 by a 
5 modulated sine wave having a second angular frequency (O 
2. With this phase modulation, the phase of the second 
exciting current is given by coot - {0 2 + m p cos(C0 2t)} 
where m p is a phase modulation index which represents 
the phase deviation amount at the maximum amplitude of 
10 the modulated wave. 

In this embodiment, b2 = b in equation (4). 
Of the magnetic field generated from the second exciting 
coil 3b when the second exciting current is supplied 
from the power supply section 4, a magnetic field 
15 component B2 which is perpendicular to both the 

electrode axis EAX and the measuring pipe axis PAX on 
the electrode axis EAX is given by 

B2 = bcos[C0 0t - {0 2 + m p cos(CD2t)}] ...(58) 
In equation ( 26 ) , bl = b2 = b, 01 = 0, and 0 
20 01 = 0 00. When the magnetic fields Bl and B2 are given 
by equations (27) and (58), we obtain 
Eac = Ec + Eve 

= bO)0rkexp{ j ( K/2 + #00)} 

+ bcoorkexpf j(-7C/2 + m p cos(C02t) 
25 + 0 2 + 0 00}] 

+ brkvexp{i( 0 00)} 

+ brkvexp[ j{m p cos( C0 2t) + 0 2 + 0 00}] 
- 40 - 
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... (59) 

The vector as the second term on the 
rig ht-hand side of equation (59). i.e.. b»0rt«Pl«-« 
,2 * mcos(»2t> * 02 ♦ 9oo>] can he rewritten to *o> 
0r,cos' (U 0t - m,cos( W 2t> - (-/2 ♦ 92 ♦ »..» as time 
expression. This time expression can aiso he rewritten 

to 

bco0rkcos{co0t - mcos(0)2t) 

_ (-7T/2 + 02 + 000)} 
= b0)0rk[cos{ coot 

(-7C/2 + 02 + 0OO)}cos{mcos(co2t)> 

+ sin{0)0t 

(-7C/2 + 02 + 0OO)}sin{mcos(co2t)}] 



(60) 



(61) 
. (62) 



5 in equation (60). cos{m p cos ( CO 2t ) } and 

sin{m p cos(C0 2t)} can be rewritten to 

P cos{m p cos(co2t)}= J 0 (m p ) + 2 £ H)"' 2 1. (m p )cos(nc»2t) 

sin{ mp cos(co2t)}= 2 £ H)<->' 2 J n (m p )cos(nco2t) 

n=l,3,... 

, 0 in equations (61) and (62), j n (m p ) (n = 0, 1. 

2 . . . . ) i. known as a Bessel function of fractional 
order. The Bessel function of fractional order. j n (m p ) 
is given by 

J n ( m p)-2- k!(n + k)!l, 2 J 
25 in equation (63), k! indicates the factorial 

of k. When only a case wherein n = 0 or 1 is applied to 



K J rm„Y +2k ...(63) 

H/n 4- 

k=0 
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equations (61) and (62), equation (60) can be rewritten 
to 

bco0rkcos{ coot - m p cos(C0 2t) 

- (-7T/2 +02+000)} 

5 = bC00rk[cos{C00t - (-7C/2 + 02 + 0OO)}J Q (m p ) 

+ sin{CO0t 

- (-7T/2 + 02 + 0 00) }2J 1 (m p )cos( 0)2t) ] 

= bco0rk[ J 0 (m p )cos{o;Ot - (- 71/2 + 0 2 + 0 00)} 
+ J^injsinf ( COO + C02)t - (-7Z/2 + 02 
10 + 000)} 

+ J 1 (m p ) sin{ ( CO 0 - C02)t - (-7t/2 + 02 
+ 000)}] 

= bCOOrkJ 0 (m p )cos{COOt - (-71/2 + 0 2 + 0 00)} 
+ bC00rkJ 1 (m p )cos{ ( C00 + C02)t - ( 02 + 000)} 
15 + bC00rkJ 1 (m p )cos{ ( COO - C02)t - ( 02 + 000)} 

. . . (64) 

As is apparent from equation (64), the second 
term on the right-hand side of equation (59) forms a 
vector bcoOrkJ 0 (m p )exp{j ( - 71/ 2 + 02 + 0 00)} on a complex 

20 plane based on the angular frequency coo. The second 

term also forms a vector bC00rkJ 1 (m p )exp{ j( 0 2 + 000)} on 
each of a complex plane based on the angular frequency 
(coo + co 2) and a complex plane based on the angular 
frequency (coo - co2). 

25 The vector as the fourth term on the 

right-hand side of equation (59), i.e., 
brkvexp [ j {rncos ( co 2t ) + 0 2 + 0 00}] can be rewritten to 
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brkvcos{C0 0t - m p cos(C0 2t) - (02 + 0 00)} as time 
expression. This time expression can also be rewritten 

brkvcos{coot - m p cos( co 2t) - (02 + 0 00)} 
5 = brkv[cos{ coot 

- (02 + 9 00) }cos{racos( C0 2t) } 

+ sin{CO0t - ( 9 2 + 9 00) }sin{m p cos( C0 2t ) } ] 

. . . (65) 

Like the second term on the right-hand side of 
10 equation (59), when the Bessel function of fractional 
order J n (ni p ) is applied, equation (65) can be rewritten 
to 

brkvcos{C0 0t - m p cos(C0 2t) - (0 2 + 0 00)} 
= brkv[cos{ 0)0t - (02 + 0OO)}J o m p 
15 + sin{C00t - (02 + 0 00) }2J 1 (m p )cos( C02t ) ] 

= brkv[ J Q (m p )cos{ coot - ( 02 + 0 00)} 

+ J 1 (m p )cos{ ( COO + C02)t - (TC/2 + 0 2 + 000)} 
+ J 1 (m p )cos{ ( C00 - C02)t - (7T/2 + 02 + 000)}] 

. • . (66) 

20 As is apparent from equation (66), the fourth 

term on the right-hand side of equation (59) forms a 
vector brkvJ 0 (m p )exp{ j ( 0 2 + 000)} on the complex plane 
based on the angular frequency coo. The fourth term 
also forms a vector brkvJ^in^expi j ( rc /2 + 0 2 + 0 00)} on 

25 each of the complex plane based on the angular frequency 
(coo + co 2) and the complex plane based on the angular 
frequency ( co 0 - co 2 ) . 
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As is apparent from the above description, the 
second and fourth terms on the right-hand side of 
equation (59) form a complex vector Epm on each of the 
complex plane based on the angular frequency (coo + co2) 
5 and the complex plane based on the angular frequency ( co 
0 - 0)2). The complex vector Epm is given by 
Epm = bO)0rkJ 1 (m p )exp{j( 0 2 + 0 00)} 

+ brkvJ 1 (m p )exp{ j( K / 2 +02+000)} 

... (67) 

10 In addition, the first to fourth terms on the 

right-hand side of equation (59) form a complex vector 
Eor on the complex plane based on the angular frequency 
coo. The complex vector Eor is given by 
Eor = bco0rkexp{ j ( K/2 + 0 00)} 
15 + brkvexp{ j( 0 00) } 

+ bC0 0rkJ Q (m p )exp{ j(-7T/2 + 0 2 + 0 00)} 
+ brkvJ Q (m p )exp{ j( 0 2 + 000)} ...(68) 
The magnitude |Epm| of the complex vector Epm 
given by equation (67) is given by 

2 2 2 1/2 

20 | Epm | = brk{J } (m p ) {COO +v)} ,..(69) 

The magnitude |Eor| of the complex vector Eorm 
given by equation (68) is given by 

| Eor | = brk[{v + J Q ( m p ) sin ( 0 2 ) CO 0 
+ J 0 (m p )cos( 0 2)v} 2 
25 + {COO - J 0 (m p )cos( 0 2) C00 

2 1/2 

+ J 0 (m p )sin( 0 2)v} ] 

. . . (70) 
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A ratio Rpm of the magnitude |Eor| of the 
complex vector Eor to the magnitude |Epm| of the complex 

vector Epm is given by 

. . i ...(71) 
Rpm = |Eor|/|Epm| 

5 The magnitude |Epm| given by equation (69) and 

the magnitude |Eor| given by equation (70) have no term 
containing the angle 0 00 . The magnitudes |Eor| and 
|Epm| contain brk. However, when the ratio Rpm is 
obtained by equation (71). brJc is erased from the ratio 

10 Rpm. 

Hence, even when the angle 0 00 changes, or 
the amplitude (the amplitude of the carrier component of 
the magnetic field B2) b of the magnetic field Bl 
shifts, the ratio Rpm does not change. When the flow 
15 rate is detected on the basis of the ratio Rpm, a flow 
rate measurement error due to the variation in delay of 
the magnetic field with respect to the exciting current 
or the shift of the amplitude of the magnetic field can 
be quickly and automatically canceled. 

When equations (69) and (70) are substituted 
into equation (71), and equation (71) is solved about v, 

we obtain 

v = [ O)O{-2J 0 (m p )sin( 6 2) 

+ {2J 0 (m p ) 2 - J 0 (m p ) 4 + 2J 0 (m p ) 2 J 1 (m p )W - 1 

2 2 4 4 -i 1/2 "» 1 

25 + 2J 1 (m p )Rpm - J x (m p ) Rpm } }) 

h.u^i + 2J.(m )cos( 6 2) - J,(m o ) 2 Rpm 2 } 

/ \ " 0 V "' p i ■ - u p 

...(72) 



20 
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(73) 



From equation (24), equation (72) can be 

rewritten to 

V= a X [ coO{-2J 0 (m p )sin( 6 2) 

+ {2J 0 (m p) 2 - J 0 (m p ) 4 + 2J 0 (m p ) 2 J 1 ("'p)W " 1 
5 + 2J 1 (m p ) Rpm - J^) Rpm } }] 

/{ j 0 (m p ) 2 + 1 + 2J o (m p )cos(0 2) - J 1 (m p )W> 

for a - 1/7 

A signal conversion section 5 
frequency- separates an interelectrode electromotive 
10 force Eac by using a filter, as in the first embodiment, 
to obtain the amplitude of the component having the 
angular frequency (coo + 0)2) or (coo - o>2> (the 
magnitude |Epm| of the complex vector Epm) . The signal 
conversion section 5 also obtains the amplitude of the 
15 component having the angular frequency coO (the 

magnitude | Bor | of the complex vector Eor) . The signal 
conversion section 5 calculates the ratio Rpm of the 

magnitudes | Eor | to | Epm | . 

A flow rate output section 6 calculates a flow 
20 velocity V of the fluid to be measured by using equation 
(73) on the basis of the ratio Rpm obtained by the 
signal conversion section 5. With the above 
arrangement, the same effect as in the first embodiment 
can be obtained. 
25 [Fifth Embodiment] 

The fifth embodiment of the present invention 

will be described next. The arrangement of an 
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electromagnetic flowmeter according to this embodiment 
is the same as in the first embodiment and will be 
described with reference to Fig. 4. A power supply 
section 4 according to this embodiment supplies a first 
5 exciting current to a first exciting coil 3a. The first 
exciting current is obtained by phase-modulating a sine 
wave carrier having a first angular frequency 0)0 by a 
modulated sine wave having a second angular frequency 0) 
2. In this embodiment, bl = b, and 6 1 = 0 in equation 
10 (3). 

With this phase modulation, the phase of the 
first exciting current is given by coot - mcos(C0 2t) + 

Of the magnetic field generated from the first 
15 exciting coil 3a when the first exciting current is 
supplied from the power supply section 4, a magnetic 
field component Bl which is perpendicular to both an 
electrode axis EAX and a measuring pipe axis PAX on the 
electrode axis EAX is given by 
20 Bl = bcos{C00t - m p cos(C02t) + it } ...(74) 

The power supply section 4 also supplies a 
second exciting current to a second exciting coil 3b. 
The second exciting current is obtained by 
phase -modulating a sine wave carrier having the same 
25 angular frequency 0)0 as that of the first exciting 

current and a predetermined phase difference 6 2 by a 
modulated sine wave having the same angular frequency 0) 
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2 as that of the modulated wave component of the first 
exciting current and an opposite phase. 

With this phase modulation, the phase of the 
second exciting current is given by coot - {9 2 + 
5 mcos(C02t)}. 

In this embodiment, b2 = b in equation (4). 
Of the magnetic field generated from the second exciting 
coil 3b when the second exciting current is supplied 
from the power supply section 4, a magnetic field 
10 component B2 which is perpendicular to both the 

electrode axis EAX and the measuring pipe axis PAX on 
the electrode axis EAX is given by equation (58). 

The magnitude |Epm| of a complex vector Epm, 
which is obtained in accordance with the same procedures 
15 as in the fourth embodiment, is given by 

1/2 2 2 2 2 

| Epm | = 2 brktJ^nip) {v + 0)0 + 0)0cos(#2) 

2 1/2 

- vcos(#2) - 2v0)0sin( 0 2) }] 

...(75) 

The magnitude | Eor | of a complex vector Eor is 

20 given by 

1/2 2 2 2 2 

|Eor| = 2 brk[J 0 (m p ) {v + to 0 - O)0cos(<9 2) 

2 1/2 

- vcos(02) + 2v0)0sin( 9 2) }] 

... (76) 

The magnitude |Epm| given by equation (75) and 
25 the magnitude |Eor| given by equation (76) have no term 
containing an angle #00. 

The magnitudes |Eor| and |Epm| contain brk. 
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However, when a ratio Rpm is obtained by equation (71), 
brk is erased from the ratio Rpm. Hence, when the flow 
rate is detected on the basis of the ratio Rpm, a flow 
rate measurement error due to the variation in delay of 
5 the magnetic field with respect to the exciting current 
or the shift of the amplitude of the magnetic field can 
be quickly and automatically canceled. 

When equations (75) and (76) are substituted 
into equation (71). and equation (71) is solved about v. 

10 we obtain 

v = coO[-{J 0 (m p ) 2 cos( 0 2)sin( 9 2) 
+ j i (m p ) 2 sin( 0 2)cos( 0 2)Rpm 2 
+ J i( m p ) 2 sin(0 2)Rpm 2 + J Q (m p ) 2 sin( 6 2 ) } 
+ 2|J,(m t )J 1 (m p ){cos(0 2) + l}Rpm| ] 
15 /{2J 0 (m p ) 2 cos( 6 2) + J 0 (m p ) 



+ j o (m p ) 2 cos( 9 2) 2 - J^m^Rpm 

2 



+ J i (m p ) 2 cos( 6 2) Rpm } 

From equation (24), equation (77) can be 

20 rewritten to 

v = a X coO[-{J 0 (m p ) 2 cos( <92)sin( 6 2) 

+ J i (m p ) 2 sin( 02)cos( 0 2)Rpm 2 

+ j i (m p ) 2 sin( 0 2)Rpm 2 + J 0 (m p ) 2 sin( 6 2) } 

+ 2 I J„(m )J,(m o ){cos( 6 2) + l}Rpm| ] 

' 0 p 1 v 

/{2J 0 (m p ) 2 cos( 9 2) + J 0 (m p ) 2 
T^\ J .ncff)?^ 2 - J (m ) 2 Rpm 2 

+ J (m ) 2 cos( 6 2) 2 Rpm 2 } 

i p 



(77) 



25 



- 49 - 



for a = 1/ r ... (78) 

A signal conversion section 5 
frequency- separates an interelectrode electromotive 
force Eac by using a filter, as in the first embodiment, 
5 to obtain the amplitude of the component having the 
angular frequency (coo + 0)2) or (coo - <x> 2 ) ( the 
magnitude |Epm| of the complex vector Epm) . 

The signal conversion section 5 also obtains 
the amplitude of the component having the angular 

10 frequency 0)0 (the magnitude |Eor| of the complex vector 
Eor) . The signal conversion section 5 calculates the 
ratio Rpm of the magnitudes |Eor| to |Epm|. 

A flow rate output section 6 calculates a flow 
velocity V of the fluid to be measured by using equation 

15 (78) on the basis of the ratio Rpm obtained by the 
signal conversion section 5. With the above 
arrangement, the same effect as in the first embodiment 
can be obtained. 
[ Sixth Embodiment ] 

20 The sixth embodiment of the present invention 

will be described next. The arrangement of an 
electromagnetic flowmeter according to this embodiment 
is the same as in the first embodiment and will be 
described with reference to Fig. 4. A power supply 

25 section 4 according to this embodiment supplies a first 
sine wave exciting current having a first angular 
frequency coo to a first exciting coil 3a. 
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In this embodiment, bl = b, and 6 1 = 0 in 
equation (3), Of the magnetic field generated from the 
first exciting coil 3a when the first exciting current 
is supplied from the power supply section 4, a magnetic 
5 field component Bl which is perpendicular to both an 

electrode axis EAX and a measuring pipe axis PAX on the 
electrode axis EAX is given by equation (27). 

The power supply section 4 also supplies a 
second exciting current to a second exciting coil 3b. 

10 The second exciting current is obtained by 

frequency-modulating a sine wave carrier having the same 
angular frequency coo as that of the first exciting 
current and a predetermined phase difference 9 2 by a 
modulated sine wave having a second angular frequency co 

15 2. 

With this phase modulation, the phase of the 
second exciting current is given by coot - {02 + 
m f sin(co2t)} where m f is a frequency modulation index 
which represents the frequency deviation amount at the 

20 maximum amplitude of the modulated wave. 

In this embodiment, b2 = b in equation (4). 
Of the magnetic field generated from the second exciting 
coil 3b when the second exciting current is supplied 
from the power supply section 4, a magnetic field 

25 component B2 which is perpendicular to both the 

electrode axis EAX and the measuring pipe axis PAX on 
the electrode axis EAX is given by 
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B2 = bcos[C00t - {02 + m £ sin(o>2t)}] ...(79) 
in equation (26),- bl = b2 = b, 0 1 = 0. and 0 
01 - 0 00. When the magnetic fields Bl and B2 are given 
by equations (27) and (79). we obtain 
Eac = EC + Eve 

= bw0rkexp{ j ( K/2 + 0 00)} 

+ bCOOrkexpt j(-7T/2 + m £ sin(C02t) 

+ 62 + 000}] 

+ brkvexp{j( 0 00) } 

+ brkvexp[j{m £ sin(to2t) + 0 2 + 0 00}] 

. . . (80) 

The vector as the second term on the 
right-hand side of equation (80), i.e., bC00rkexp[j{- * 
12 * m £ sin(C02t) + 0 2 + 0 00}] can be rewritten to bco 

i TT/9+ 02 + 0 00)} as time 
0rkcos{co0t - m £ sxn(co2t) - {-n/t + o z ^ , j 

expression. This time expression can also be rewritten 
to 

b<x>0rkcos{co0t - m £ sin(o>2t) 

_ (-7T/2 + 0 2 + 000)} 
= b co0rk[ cos { coot 

- (-7E/2 + 02 + 0OO)}cos{m £ sin(CO2t)} 

+ sin{CO0t 

- (-7T/2 + 02 + 0 00) }sin{m £ sin( C02t)}] 



(81) 



In equation (81), cos{m f sin( CO 2t ) } and 
sin{m f sin(0)2t)} can be rewritten to 
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cos{m f sin(co2t)} = J 0 (m f ) + 2 £<-!)» ' 2 J n (m f )sin(nc»2t) • • • (82) 

n=2,4,- 

si„{m f sin(co2t)} = 2 ]T (-l) (n_1) /2 J „ (m f ) sin(n©2t) • • • (83) 

n=l,3,... 

In equations (82) and (83), the Bessel 
function of fractional order, j n (m £ ) (n = 0, 1, 2,...) is 
given by 

T , > y (-D k (VY** •••< 84 > 

J n( m f ) = 2, k!(n + k ),X 2 J 

When only a case wherein n = 0 or 1 is applied 
to equations (82) and (83), equation (81) can be 

rewritten to 

b0>0rkcos{ wot - m £ sin(co2t) 

- (-7T/2 + 02 + 000)} 

= bW0rk[cos{C00t - (-7C/2 + 02 + 0OO)}J o (m £ ) 

+ sin{coot 

- (-7T/2 + 02 + 0OO)}2J 1 (m £ )sin((o2t)] 

= bWOrk[J 0 (m £ )cos{o>Ot - (-7C/2 + 02 + 000)} 

- J 1 (m £ )cos{(C0 0 + C02)t - (-7T/2 + 92 
+ 000)} 

+ J 1 (m £ )cos{(0)0 - 0)2)t - (-7C/2 + 02 
+ 000)}] 

. . . (85) 

As is apparent from equation (85), the second 
term on the right-hand side of equation (80) forms a 
vector bcoOrkJ 0 (m £ )exp{j(-7t/2 + 0 2 + 000)} on a complex 
plane based on the angular frequency coo, a vector bco 
0rk{-3 1 (m £ )}exp{j(-7r/2 + 0 2 + 0 00)} on a complex plane 



o 



based on the angular frequency (coo + C02), and a vector 
bC00rkj 1 (m f )exp{j(-7r/2 + 0 2 + 0 00)} on a complex plane 
based on the angular frequency (coO - co2). 

The vector as the fourth term on the 
right-hand side of equation (80), i.e., 

brkvexpt j{m f sin( C02t ) + 0 2 + 9 00}] can be rewritten to 
brkvcos{C00t - m f sin(a>2t) - ( 0 2 + 0 00)} as time 
expression. This time expression can also be rewritten 
to 

brkvcos{coot - m f sin(C0 2t) - (02 + 0 00)} 
= brkv [ cos { coot 

- (02 + 0 00) }cos{m f sin( C0 2t) } 

+ sin{co0t - ( 0 2 + 0 00) }sin{m f sin( C0 2t ) } ] 

. . . (86) 

Like the second term on the right-hand side of 
equation (80), when the Bessel function of fractional 
order J n (m f ) is applied, equation (86) can be rewritten 
to 

brkvcos{ coot - m f sin(C0 2t) - ( 0 2 + 0 00)} 
= brkv[cos{ coot - (02 + 0OO)}J o m f 

+ sin{coot - ( 0 2 + 0 00)}2J 1 (m f )sin( C02t) ] 
= brkv[ J 0 (m f )cos{co0t - ( 02 + 0 00)} 

- J 1 (m f )cos{ ( C00 + C02)t - ( 02 + 000)} 
+ J 1 (m f )cos{( 600 - C02)t - ( 02 + 000)}] 

. . . (87) 

As is apparent from equation (87), the fourth 
term on the right-hand side of equation (80) forms a 
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vector brkvJ Q (m f )exp{ j ( 9 2 + 000)} on the complex plane 
based on the angular frequency 0)0, a vector 
brkv{- j (m ) }exp{ j ( 0 2 + 000)} on the complex plane based 
on the angular frequency (coo + co2), and a vector 
brkvj 1 (m f )exp{ j( 6 2 + 000)} on the complex plane based on 
the angular frequency (too - (D2). 

As is apparent from the above description, the 
second and fourth terms on the right-hand side of 
equation (80) form complex vectors having the same 
magnitude and reverse directions on the complex plane 
based on the angular frequency (coo + C0 2) and the 
complex plane based on the angular frequency (coo - co 
2). 

A complex vector Efm formed on the complex 
plane based on the angular frequency (o>0 - 0)2) is 
given by 

Efm = bO)0rkJ 1 (m f )exp{ j(-7T/2 + 0 2 + 0 00)} 
+ brkvJ 1 (m f )exp{ j( 0 2 + 000)} 

... (88) 

In addition, the first to fourth terms on the 
right-hand side of equation (80) form a complex vector 
Eor on the complex plane based on the angular frequency 
coo. The complex vector Eor is given by 
Eor = bC0 0rkexp{ j ( rt/2 + 0 00)} 
+ brkvexp{ j ( 0 00) } 

+ bcoOrkJ 0 (m f )exp{j(-7r/2 + 02 + 000)} 

+ brkvJ 0 (m f )exp{ j( 0 2 + 0 00)} ...(89) 
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A ratio Rfm of the magnitude | Eor | of the 
complex vector Eor to the magnitude |Efm| of the complex 
vector Efm is given by 

Rfm = |Eor|/|Efm| ...(90) 
When the magnitude | Eor | of the complex vector 
Eor obtained from equation (89) and the magnitude |Efm| 
of the complex vector Efm obtained from equation (88) 
are substituted into equation (90), and equation (90) is 
solved about v, we obtain 

v = [ coO{-2J 0 (m f )sin( 9 2) 

2 2 2 

2 t / \ _i_ it (m \ .1 fm } Rfm - 1 



+ {2J 0 (m £ ) 2 - J 0 (m £ ) + 2J 0 (m £ ) J 1 (m £ ) Rfm 
+ 2J £ (m £ ) Rfm - Jj(m £ ) Rfm } }] 
/{j 0 (m £ ) 2 + 1 + 2J o (m £ )cos(0 2) - J 1 (m £ )W} 



. . . (91) 

Equation (91) has no term containing an angle 
0 00 or the amplitude (the amplitude of the carrier 
component of the magnetic field B2) b of the magnetic 
field Bl. Hence, when the flow rate is detected on the 
basis of the ratio Rfm, a flow rate measurement error 
due to the variation in delay of the magnetic field with 
respect to the exciting current or the shift of the 
amplitude of the magnetic field can be quickly and 
automatically canceled. 

From equation (24). equation (91) can be 

rewritten to 

Y= a X [ coo{-2J 0 (m £ )sin( 9 2) 

+ {2J 0 (m £ ) 2 - J 0 (m £ ) 4 + 2J 0 (m £ ) 2 J 1 (m £ )W - 1 
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2 2 4 4 1/2 

+ 2 J 1 (m f ) Rfm - Jj(m f ) Rfm } }] 
/{J 0 (m f ) 2 + 1 + 2J 0 (m f )cos( 9 2) - J 1 (m £ ) 2 Rfm 2 } 
for a = 1/7 ... (92) 

A signal conversion section 5 
5 frequency- separates an interelectrode electromotive 

force Eac by using a filter, as in the first embodiment, 
to obtain the amplitude of the component having the 
angular frequency (coo - a)2) (the magnitude |Efm| of 
the complex vector Efm) . The signal conversion section 
10 5 also obtains the amplitude of the component having the 
angular frequency coo (the magnitude |Eor| of the 
complex vector Eor) . The signal conversion section 5 
calculates the ratio Rfm of the magnitudes |Eorj to 
| Efm | . 

15 A flow rate output section 6 calculates a flow 

velocity V of the fluid to be measured by using equation 
(92) on the basis of the ratio Rfm obtained by the 
signal conversion section 5. With the above 
arrangement, the same effect as in the first embodiment 

20 can be obtained. 

[Seventh Embodiment] 

The seventh embodiment of the present 
invention will be described next. The arrangement of an 
electromagnetic flowmeter according to this embodiment 

25 is the same as in the first embodiment and will be 
described with reference to Fig. 4. 

A power supply section 4 according to this 
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embodiment supplies a first exciting current to a first 
exciting coil 3a. The first exciting current is 
obtained by frequency-modulating a sine wave carrier 
having a first angular frequency coo by a modulated sine 
5 wave having a second angular frequency C0 2. In this 
embodiment, bl = b, and 01 = 0 in equation (3). 

With this frequency modulation, the phase of 
the first exciting current is given by coot - m f sin(co2t 

+ 7t ) . 

10 Of the magnetic field generated from the first 

exciting coil 3a when the first exciting current is 
supplied from the power supply section 4, a magnetic 
field component Bl which is perpendicular to both an 
electrode axis EAX and a measuring pipe axis PAX on the 

15 electrode axis EAX is given by 

Bl = bcos{0)0t - m f sin(0)2t + 7t) } ...(93) 
The power supply section 4 also supplies a 
second exciting current to a second exciting coil 3b. 
The second exciting current is obtained by 

20 frequency-modulating a sine wave carrier having the same 
angular frequency coo as that of the first exciting 
current and a predetermined phase difference 9 2 by a 
modulated sine wave having the same angular frequency co 
2 as that of the modulated wave component of the first 

25 exciting current and an opposite phase. 

With this frequency modulation, the phase of 
the second exciting current is given by coot - 9 2 + 
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m f sin( co 2t) . 

In this embodiment, b2 = b in equation (4). 
Of the magnetic field generated from the second exciting 
coil 3b when the second exciting current is supplied 
5 from the power supply section 4, a magnetic field 
component B2 which is perpendicular to both the 
electrode axis EAX and the measuring pipe axis PAX on 
the electrode axis EAX is given by equation (79), 

The magnitude |Eor| of a complex vector Eor 
10 and the magnitude |Efm| of a complex vector Efm are 

obtained in accordance with the same procedures as in 
the sixth embodiment. When the magnitudes |Eor| and 
| Efm | are substituted into equation (90), and equation 
(90) is solved about v, we obtain 
15 v = 0)0[-{ J 0 (m f ) 2 cos( 9 2)sin( 9 2) 

+ J^m^'sinf 9 2)cos( 6 2)Rfm 2 
+ J 1 (m f ) 2 s-in( 9 2)R£m + J fl (m f ) 2 sin( 9 2) } 
+ 2| J Q (m f ) J 1 (m f ){cos( 9 2) + l}Rfm| ] 
/{2J 0 (m f ) 2 cos( 9 2) + J 0 (m f ) 2 
20 + J Q (m f ) 2 cos( 9 2) 2 - J 1 (m f ) 2 Rfm 2 

+ J 1 (m f ) 2 cos( 9 2) 2 Rfm 2 } 

. . .(94) 

Equation (94) has no term containing an angle 
9 00 or an amplitude b of the carrier components of the 
25 magnetic fields Bl and B2 . Hence, when the flow rate is 
detected on the basis of a ratio Rfm, a flow rate 
measurement error due to the variation in delay of the 
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magnetic field with respect to the exciting current or 
the shift of the amplitude of the magnetic field can be 
quickly and automatically canceled. 

From equation (24), equation (94) can be 
rewritten to 

V= a X coO[-{J 0 (m f ) 2 cos( 02)sin( 9 2) 

+ J ] (m f ) 2 sin( 9 2)cos( 9 2)Rfm 2 

+ J 1 (m f ) 2 sin( 9 2)Rfm 2 + J Q (m f ) 2 sin( 9 2 ) } 

+ 2| J Q (m f ) J 1 (m f ){cos( 9 2) + l}Rfm| ] 

/{2J 0 (m f )cos( 9 2) + J 0 (m f ) 2 

2 2 2 2 

+ J 0 (m f ) cos( 9 2) - Jjtm^ Rfm 

+ J 1 (m f ) 2 cos( 9 2) 2 Rfm 2 } 
for a = 1/ r ... (95) 

A signal conversion section 5 
frequency- separates an interelectrode electromotive 
force Eac by using a filter, as in the first embodiment, 
to obtain the amplitude of the component having the 
angular frequency (coo - CO 2) (the magnitude |Efm| of 
the complex vector Efm) . The signal conversion section 
5 also obtains the amplitude of the component having the 
angular frequency o>0 (the magnitude |Eor| of the 
complex vector Eor) . The signal conversion section 5 
calculates the ratio Rfm of the magnitudes |Eor| to 
| Efm | • 

A flow rate output section 6 calculates a flow 
velocity V of the fluid to be measured by using equation 
(95) on the basis of the ratio Rfm obtained by the 
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signal conversion section 5. With the above 
arrangement, the same effect as in the first embodiment 
can be obtained. 

In the fourth to seventh embodiments, only a 
5 case wherein n = 0 or 1 is applied in expanding the 
Bessel function of fractional order. Instead, the 
amplitude ratio Rpm or Rfm is obtained by separating a 
component having a sum frequency (coo + |cl>2) or a 
difference frequency (0)0 - 6 C0 2) of the first angular 

10 frequency 0)0 and a third angular frequency (an angular 
frequency £ -times higher than the second angular 
frequency CD 2) |co2 from the interelectrode 
electromotive force Eac. The flow velocity V of the 
fluid to be measured is calculated from the amplitude 

15 ratio Rpm or Rfm. In this case, even when 4 is an 
integer of 2 or more, the flow velocity V can be 
calculated, as in the fourth to seventh embodiments, by 
applying a case wherein n = 2 or more in expanding the 
Bessel function of fractional order. 

20 In the first to seventh embodiments, in-phase 

component noise can be removed. Hence, the rectangular 
wave exciting method need not be used. Since the sine 
wave exciting method which uses a sine wave for an 
exciting current can be used, high-frequency excitation 

25 can be executed. When high-frequency excitation is 

used, 1/f noise can be removed, and the response to a 
change in flow rate can be increased. 
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The magnetic fields Bl and B2 applied to the 
fluid to be measured only need to satisfy the conditions 
described in each embodiment. Hence, the exciting coils 
may be arranged symmetrically by setting them such that 
5 the offset distance dl from the plane PLN to the axis of 
the first exciting coil 3a is equal to the offset 
distance 6.2 for the plane PLN to the axis of the second 
exciting coil 3b. Alternatively, the exciting coils may 
be arranged asymmetrically by setting different offset 

10 distances dl and d2. 

The angle made by the electrode axis EAX and 
the axis of the first exciting coil 3a may be equal to 
the angle made by the electrode axis EAX and the axis of 
the second exciting coil 3b (e.g., 90°). Alternatively, 

15 as shown in Fig. 8, the angle made by the electrode axis 
EAX and the axis of the first exciting coil 3a may be 
different from the angle made by the electrode axis EAX 
and the axis of the second exciting coil 3b. In the 
first to seventh embodiments, the phase difference 9 2 

20 may be 0. 

As the electrodes 2a and 2b used in the first 
to seventh embodiments, electrodes which are exposed 
from the inner wall of the measuring pipe 1 and come 
into contact with the fluid to be measured may be used, 
25 as shown in Fig. 9. Alternatively, as shown in Fig. 10, 
capacitively coupled electrodes which do not come into 
contact with the fluid to be measured may be used. When 
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capacitively coupled electrodes are used, the electrodes 
2a and 2b are covered with a lining 10 made of ceramic 
or Teflon (registered trademark) and formed on the inner 
wall of the measuring pipe 1 . 
5 In the first to seventh embodiments, the two 

electrodes 2a and 2b are used. However, the present 
invention is not limited to this. The present invention 
can also be applied to an electromagnetic flowmeter 
having only one electrode. When only one electrode is 

10 used, a ground ring is provided in the measuring pipe 1 
to set the potential of the fluid to be measured to the 
ground potential. An electromotive force (the potential 
difference from the ground potential) generated in the 
single electrode is detected by the signal conversion 

15 section 5, When the two electrodes 2a and 2b are used, 
the electrode axis EAX forms a straight line that 
connects the electrodes 2a and 2b. When only one 
electrode is used, it is assumed that a virtual 
electrode is arranged on the plane PLN including the 

20 single real electrode at a position opposite to the real 
electrode with respect to the measuring pipe axis PAX. 
The straight line that connects the real electrode and 
the virtual electrode at this time corresponds to the 
electrode axis EAX. 

25 The means for calculating the ratio Ram in the 

signal conversion section 5 and the flow rate output 
section 6 in the first embodiment, the means for 
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obtaining the phase differences <£>or and 4) am in the 
signal conversion section 5 and the flow rate output 
section 6 in the second embodiment, the means for 
calculating the ratio Ram in the signal conversion 
5 section 5 and the flow rate output section 6 in the 
third embodiment, the means for calculating the ratio 
Rpm in the signal conversion section 5 and the flow rate 
output section 6 in the fourth and fifth embodiments, 
and the means for calculating the ratio Rfm in the 

10 signal conversion section 5 and the flow rate output 
section 6 in the sixth and seventh embodiments can be 
implemented by, e.g., a computer. 

According to the present invention, a first 
magnetic field having a first frequency is applied to a 

15 fluid while a second magnetic field obtained by 

amplitude-, phase-, or frequency-modulating a carrier 
having the first frequency by a modulated wave having a 
second frequency is simultaneously applied to the fluid. 
Alternatively, a first magnetic field obtained by 

20 amplitude-, phase-, or frequency-modulating a carrier 
having a first frequency by a modulated wave having a 
second frequency is applied to a fluid while a second 
magnetic field obtained by amplitude-, phase-, or 
frequency-modulating a carrier having the first 

2 5 frequency by a modulated wave having the same frequency 
as that of the modulated wave of the first magnetic 
field and an opposite phase is simultaneously applied to 
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the fluid. 

Accordingly, a plurality of frequency 
components, i.e., the first frequency and the sum and 
difference frequencies of the first frequency and a 
5 third frequency (a frequency as an integer multiple of 
the second frequency) are generated in the 
interelectrode electromotive force. From two of the 
plurality of frequency components, an asymmetrical 
exciting characteristic parameter (amplitude ratio or 

10 phase difference) can be obtained, which depends on the 
flow rate of the fluid and does not depend on the 
variation in delay of the magnetic field with respect to 
the exciting current or the shift of the amplitude of 
the magnetic field. 

15 When this asymmetrical exciting characteristic 

parameter is used, a flow rate can be calculated while 
quickly and automatically canceling a flow rate 
measurement error due to the variation in delay of the 
magnetic field with respect to the exciting current or 

20 the shift of the amplitude of the magnetic field. 

In addition, since in- phase component noise 
can be removed, the rectangular wave exciting method 
need not be used. Since the sine wave exciting method 
can be used, high-frequency excitation can be executed. 

25 As a result, the flow rate can be accurately measured. 
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